Exome-wide association study to identify rare variants influencing COVID-19 outcomes: Results from the Host Genetics Initiative.

Host genetics is a key determinant of COVID-19 outcomes. Previously, the COVID-19 Host Genetics Initiative genome-wide association study used common variants to identify multiple loci associated with COVID-19 outcomes. However, variants with the largest impact on COVID-19 outcomes are expected to be rare in the population. Hence, studying rare variants may provide additional insights into disease susceptibility and pathogenesis, thereby informing therapeutics development. Here, we combined whole-exome and whole-genome sequencing from 21 cohorts across 12 countries and performed rare variant exome-wide burden analyses for COVID-19 outcomes. In an analysis of 5,085 severe disease cases and 571,737 controls, we observed that carrying a rare deleterious variant in the SARS-CoV-2 sensor toll-like receptor TLR7 (on chromosome X) was associated with a 5.3-fold increase in severe disease (95% CI: 2.75-10.05, p = 5.41x10-7). This association was consistent across sexes. These results further support TLR7 as a genetic determinant of severe disease and suggest that larger studies on rare variants influencing COVID-19 outcomes could provide additional insights.

A Timeframe for SARS-CoV-2 Genomes: A Proof of Concept for Postmortem Interval Estimations.

Establishing the timeframe when a particular virus was circulating in a population could be useful in several areas of biomedical research, including microbiology and legal medicine. Using simulations, we demonstrate that the circulation timeframe of an unknown SARS-CoV-2 genome in a population (hereafter, estimated time of a queried genome [QG]; tE-QG) can be easily predicted using a phylogenetic model based on a robust reference genome database of the virus, and information on their sampling dates. We evaluate several phylogeny-based approaches, including modeling evolutionary (substitution) rates of the SARS-CoV-2 genome (~10-3 substitutions/nucleotide/year) and the mutational (substitutions) differences separating the QGs from the reference genomes (RGs) in the database. Owing to the mutational characteristics of the virus, the present Viral Molecular Clock Dating (VMCD) method covers timeframes going backwards from about a month in the past. The method has very low errors associated to the tE-QG estimates and narrow intervals of tE-QG, both ranging from a few days to a few weeks regardless of the mathematical model used. The SARS-CoV-2 model represents a proof of concept that can be extrapolated to any other microorganism, provided that a robust genome sequence database is available. Besides obvious applications in epidemiology and microbiology investigations, there are several contexts in forensic casework where estimating tE-QG could be useful, including estimation of the postmortem intervals (PMI) and the dating of samples stored in hospital settings.

Effectiveness of COVID-19 vaccine booster in the general population and in subjects with comorbidities. A population-based study in Spain.

BACKGROUND: Research on the effectiveness of COVID-19 booster-based vaccine schedule is ongoing and real-world data on vaccine effectiveness (VE) in comorbid patients are limited. We aimed to estimate booster dose VE against SARS-CoV-2 infection and COVID-19 severity in the general population and in comorbid patients. METHOD: A retrospective test-negative control study was undertaken in Galicia-Spain (December 2020-November 2021). VE and 95% confidence interval (CI) were estimated using multivariate logistic regression models. RESULTS: 1,512,415 (94.13%) negative and 94,334 (5.87%) positive SARS-CoV-2 test results were included. A booster dose of COVID-19 vaccine is associated with substantially higher protection against SARS-CoV-2 infection than vaccination without a booster [VEboosted = 87% (95%CI: 83%; 89%); VEnon-boosted = 66% (95%CI: 65%; 67%)]. The high VE was observed in all ages, but was more pronounced in subjects older than 65 years. VE against COVID-19 severity was analyzed in a mixed population of boosted and non-boosted individuals and considerable protection was obtained [VE: hospitalization = 72% (95%CI: 68%; 75%); intensive care unit administration = 83% (95%CI: 78%; 88%), in-hospital mortality = 66% (95%CI: 53%; 75%)]. Boosted comorbid patients are more protected against SARS-CoV-2 infection than those who were non-boosted. This was observed in a wide range of major diseases including cancer (81% versus 54%), chronic obstructive pulmonary disease (84% versus 61%), diabetes (84% versus 65%), hypertension (82% versus 65%) and obesity (91% versus 67%), among others. CONCLUSIONS: A booster dose of COVID-19 vaccine increases the protection against SARS-CoV-2 infection and COVID-19 severity in the general population and in comorbid patients.

Evaluation of BNT162b2 Vaccine Effectiveness in Galicia, Northwest Spain.

Investigating vaccine effectiveness (VE) in real-world conditions is crucial, especially its variation across different settings and populations. We undertook a test-negative control study in Galicia (Northwest Spain) to assess BNT162b2 effectiveness against acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection as well as COVID-19 associated hospitalization, intensive care unit (ICU) admission and mortality. A total of 44,401 positive and 817,025 negative SARS-CoV-2 test results belonging to adults were included. Adjusted odds ratios of vaccination and their 95% confidence interval (CI) were estimated using multivariate logistic-regression models. BNT162b2 showed high effectiveness in reducing SARS-CoV-2 infections in all age categories, reaching maximum VE ≥ 14 days after administering the second dose [18-64 years: VE = 92.9% (95%CI: 90.2-95.1); 65-79 years: VE = 85.8% (95%CI: 77.3-91.9), and ≥80 years: VE = 91.4% (95%CI: 87.9-94.1)]. BNT162b2 also demonstrated effectiveness in preventing COVID-19 hospitalization for all age categories, with VE more pronounced for those aged ≥80 years [VE = 60.0% (95%CI: 49.4-68.3)]. Moreover, there was a considerable reduction in ICU admission [VE = 88.0% (95%CI: 74.6-95.8)] and mortality [VE = 38.0% (95%CI: 15.9-55.4)] in the overall population. BNT162b2 showed substantial protection against SARS-CoV-2 infections and COVID-19 severity. Our findings would prove useful for systematic reviews and meta-analysis on COVID-19 VE.

Evaluation of BNT162b2 Vaccine Effectiveness in Galicia, Northwest Spain

Investigating vaccine effectiveness (VE) in real-world conditions is crucial, especially its variation across different settings and populations. We undertook a test-negative control study in Galicia (Northwest Spain) to assess BNT162b2 effectiveness against acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection as well as COVID-19 associated hospitalization, intensive care unit (ICU) admission and mortality. A total of 44,401 positive and 817,025 negative SARS-CoV-2 test results belonging to adults were included. Adjusted odds ratios of vaccination and their 95% confidence interval (CI) were estimated using multivariate logistic-regression models. BNT162b2 showed high effectiveness in reducing SARS-CoV-2 infections in all age categories, reaching maximum VE ≥14 days after administering the second dose [18–64 years: VE = 92.9% (95%CI: 90.2–95.1);65–79 years: VE = 85.8% (95%CI: 77.3–91.9), and ≥80 years: VE = 91.4% (95%CI: 87.9–94.1)]. BNT162b2 also demonstrated effectiveness in preventing COVID-19 hospitalization for all age categories, with VE more pronounced for those aged ≥80 years [VE = 60.0% (95%CI: 49.4–68.3)]. Moreover, there was a considerable reduction in ICU admission [VE = 88.0% (95%CI: 74.6–95.8)] and mortality [VE = 38.0% (95%CI: 15.9–55.4)] in the overall population. BNT162b2 showed substantial protection against SARS-CoV-2 infections and COVID-19 severity. Our findings would prove useful for systematic reviews and meta-analysis on COVID-19 VE.

A multi-tissue study of immune gene expression profiling highlights the key role of the nasal epithelium in COVID-19 severity.

Coronavirus Disease-19 (COVID-19) symptoms range from mild to severe illness; the cause for this differential response to infection remains unknown. Unravelling the immune mechanisms acting at different levels of the colonization process might be key to understand these differences. We carried out a multi-tissue (nasal, buccal and blood; n = 156) gene expression analysis of immune-related genes from patients affected by different COVID-19 severities, and healthy controls through the nCounter technology. Mild and asymptomatic cases showed a powerful innate antiviral response in nasal epithelium, characterized by activation of interferon (IFN) pathway and downstream cascades, successfully controlling the infection at local level. In contrast, weak macrophage/monocyte driven innate antiviral response and lack of IFN signalling activity were present in severe cases. Consequently, oral mucosa from severe patients showed signals of viral activity, cell arresting and viral dissemination to the lower respiratory tract, which ultimately could explain the exacerbated innate immune response and impaired adaptative immune responses observed at systemic level. Results from saliva transcriptome suggest that the buccal cavity might play a key role in SARS-CoV-2 infection and dissemination in patients with worse prognosis. Co-expression network analysis adds further support to these findings, by detecting modules specifically correlated with severity involved in the abovementioned biological routes; this analysis also provides new candidate genes that might be tested as biomarkers in future studies. We also found tissue specific severity-related signatures mainly represented by genes involved in the innate immune system and cytokine/chemokine signalling. Local immune response could be key to determine the course of the systemic response and thus COVID-19 severity. Our findings provide a framework to investigate severity host gene biomarkers and pathways that might be relevant to diagnosis, prognosis, and therapy.

Superspreading in the emergence of COVID-19 variants.

Superspreading and variants of concern (VOC) of the human pathogen severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are the main catalyzers of the coronavirus disease 2019 (COVID-19) pandemic. However, measuring their individual impact is challenging. By examining the largest database of SARS-CoV-2 genomes The Global Initiative on Sharing Avian Influenza Data [GISAID; n >1.2 million high-quality (HQ) sequences], we present evidence suggesting that superspreading has had a key role in the epidemiological predominance of VOC. There are clear signatures in the database compatible with large superspreading events (SSEs) coinciding chronologically with the worst epidemiological scenarios triggered by VOC. The data suggest that, without the randomness effect of the genetic drift facilitated by superspreading, new VOC of SARS-CoV-2 would have had more limited chance of success.

Seroprevalence of SARS-CoV-2 Among Pediatric Healthcare Workers in Spain.

Spain is one of the countries most severely affected by the SARS-CoV-2 pandemic, with almost 190,000 cases as of April 18, 2020. As healthcare workers (HCW) are one of the groups hardest hit by the infection, it is important to know the seroprevalence of antibodies against SARS-CoV-2 in pediatric departments. We performed 175 immunoglobulin (Ig)M and IgG immunochromatographic rapid tests in the personnel working at the Pediatric Department of the Hospital Clínico Universitario of Santiago de Compostela (Spain), including pediatricians, residents, nurses, and other staff, on days 31-33 since the lockdown started. Seven out of the 175 tests were positive, including four for IgM and three for IgG, leading to a seroprevalence of 4.0% (95% CI: 1.1-6.9%). Only one of them had symptoms at the time of testing (sore throat). All seropositive cases yielded negative RT-PCR of the upper and lower respiratory tract. This is the first SARS-CoV-2 serological survey among HCWs reported in Spain. Notwithstanding the test limitations, our results reveal that personal protection policy and lockdown measures have been effective to limit population exposure. The low seroprevalence rate poses a significant challenge for the next strategic steps of pandemic control.

CovidPhy: A tool for phylogeographic analysis of SARS-CoV-2 variation.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the pathogen responsible for the coronavirus disease 2019 (COVID-19) pandemic. SARS-CoV-2 genomes have been sequenced massively and worldwide and are now available in different public genome repositories. There is much interest in generating bioinformatic tools capable to analyze and interpret SARS-CoV-2 variation. We have designed CovidPhy (http://covidphy.eu), a web interface that can process SARS-CoV-2 genome sequences in plain fasta text format or provided through identity codes from the Global Initiative on Sharing Avian Influenza Data (GISAID) or GenBank. CovidPhy aggregates information available on the large GISAID database (>1.49 M genomes). Sequences are first aligned against the reference sequence and the interface provides different sources of information, including automatic classification of genomes into a pre-computed phylogeny and phylogeographic information, haplogroup/lineage frequencies, and sequencing variation, indicating also if the genome contains known variants of concern (VOC). Additionally, CovidPhy allows searching for variants and haplotypes introduced by the user and includes a list of genomes that are good candidates for being responsible for large outbreaks worldwide, most likely mediated by important superspreading events, indicating their possible geographic epicenters and their relative impact as recorded in the GISAID database.

Pitfalls of barcodes in the study of worldwide SARS-CoV-2 variation and phylodynamics.

Analysis of SARS-CoV-2 genome variation using a minimal number of selected informative sites conforming a genetic barcode presents several drawbacks. We show that purely mathematical procedures for site selection should be supervised by known phylogeny (i) to ensure that solid tree branches are represented instead of mutational hotspots with poor phylogeographic proprieties, and (ii) to avoid phylogenetic redundancy. We propose a procedure that prevents information redundancy in site selection by considering the cumulative informativeness of previously selected sites (as a proxy for phylogenetic-based criteria). This procedure demonstrates that, for short barcodes (e.g., 11 sites), there are thousands of informative site combinations that improve previous proposals. We also show that barcodes based on worldwide databases inevitably prioritize variants located at the basal nodes of the phylogeny, such that most representative genomes in these ancestral nodes are no longer in circulation. Consequently, coronavirus phylodynamics cannot be properly captured by universal genomic barcodes because most SARS-CoV-2 variation is generated in geographically restricted areas by the continuous introduction of domestic variants.

Mapping genome variation of SARS-CoV-2 worldwide highlights the impact of COVID-19 super-spreaders.

The human pathogen severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the major pandemic of the twenty-first century. We analyzed more than 4700 SARS-CoV-2 genomes and associated metadata retrieved from public repositories. SARS-CoV-2 sequences have a high sequence identity (>99.9%), which drops to >96% when compared to bat coronavirus genome. We built a mutation-annotated reference SARS-CoV-2 phylogeny with two main macro-haplogroups, A and B, both of Asian origin, and more than 160 sub-branches representing virus strains of variable geographical origins worldwide, revealing a rather uniform mutation occurrence along branches that could have implications for diagnostics and the design of future vaccines. Identification of the root of SARS-CoV-2 genomes is not without problems, owing to conflicting interpretations derived from either using the bat coronavirus genomes as an outgroup or relying on the sampling chronology of the SARS-CoV-2 genomes and TMRCA estimates; however, the overall scenario favors haplogroup A as the ancestral node. Phylogenetic analysis indicates a TMRCA for SARS-CoV-2 genomes dating to November 12, 2019, thus matching epidemiological records. Sub-haplogroup A2 most likely originated in Europe from an Asian ancestor and gave rise to subclade A2a, which represents the major non-Asian outbreak, especially in Africa and Europe. Multiple founder effect episodes, most likely associated with super-spreader hosts, might explain COVID-19 pandemic to a large extent.

Phylogeography of SARS-CoV-2 pandemic in Spain: a story of multiple introductions, micro-geographic stratification, founder effects, and super-spreaders.

Spain has been one of the main global pandemic epicenters for coronavirus disease 2019 (COVID-19). Here, we analyzed >41 000 genomes (including >26 000 high-quality (HQ) genomes) downloaded from the GISAID repository, including 1 245 (922 HQ) sampled in Spain. The aim of this study was to investigate genome variation of novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and reconstruct phylogeographic and transmission patterns in Spain. Phylogeographic analysis suggested at least 34 independent introductions of SARS-CoV-2 to Spain at the beginning of the outbreak. Six lineages spread very successfully in the country, probably favored by super-spreaders, namely, A2a4 (7.8%), A2a5 (38.4%), A2a10 (2.8%), B3a (30.1%), and B9 (8.7%), which accounted for 87.9% of all genomes in the Spanish database. One distinct feature of the Spanish SARS-CoV-2 genomes was the higher frequency of B lineages (39.3%, mainly B3a+B9) than found in any other European country. While B3a, B9, (and an important sub-lineage of A2a5, namely, A2a5c) most likely originated in Spain, the other three haplogroups were imported from other European locations. The B3a strain may have originated in the Basque Country from a B3 ancestor of uncertain geographic origin, whereas B9 likely emerged in Madrid. The time of the most recent common ancestor (TMRCA) of SARS-CoV-2 suggested that the first coronavirus entered the country around 11 February 2020, as estimated from the TMRCA of B3a, the first lineage detected in the country. Moreover, earlier claims that the D614G mutation is associated to higher transmissibility is not consistent with the very high prevalence of COVID-19 in Spain when compared to other countries with lower disease incidence but much higher frequency of this mutation (56.4% in Spain vs. 82.4% in rest of Europe). Instead, the data support a major role of genetic drift in modeling the micro-geographic stratification of virus strains across the country as well as the role of SARS-CoV-2 super-spreaders.